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An exact diagonalization study reveals that a matter-wave bright soliton and the Goldstone 
mode are simultaneously created in a quasi-one-dimensional attractive Bose-Einstein condensate by 
superpositions of quasi-degenerate low-lying many-body states. Upon formation of the soliton the 
maximum eigenvalue of the single-particle density matrix increases dramatically, indicating that 
a fragmented condensate converts into a single condensate as a consequence of the breaking of 
translation symmetry. 

PACS numbers: 03.75.Hh,03.75.Lm 



Fragmentation of a Bose-Einstcin condensate (BEC), 
which occurs as a consequence of a certain exact symme- 
try of the systemj^has recently been discussed in a number 
of articles 0,Qjy|- ^^ contrast to the conventional BEC, 
characterized by a unique macroscopic eigenvalue in the 
single-particle density matrix [j|, the fragmented BEC 
is characterized by more than one macroscopic eigen- 
value m . If the system has an exact symmetry and if the 
many-body theory predicts fragmentation of the ground 
state, the Gross-Pitaevskii (GP) mean-field theory does 
not predict a fragmented condensate but approximates 
it with a single condensate whose symmetry is sponta- 
neously broken. For example, a quasi-one-dimensional 
(ID) BEC with attractive interaction forms bright soli- 
tons [3, which are well described by the GP theory [y|. 
Efforts to elucidate how such symmetry-broken states 
emerge from exact many-body states have been made 
in diverse systems [7|. 

In this Letter, we show that the formation of a broken- 
symmetry soliton and an enhancement of the condensate 
fraction arc caused by superpositions of the low-lying 
states of the symmetry-preserving many-body Hamilto- 
nian. We find that the many-body spectrum exhibits a 
number of quasi-degenerate states in the regime where 
the exact ground state is a fragmented condensate. Su- 
perposition of these quasi-degenerate levels simultane- 
ously generates the broken-symmetry bright soliton and 
the Goldstone mode, accompanied by a significant in- 
crease in the condensate fraction. By introducing a small 
symmetry-breaking perturbation or by considering the 
action of a quantum measurement, we explicitly show 
that the fragmented condensate is very fragile against the 
soliton formation. Also elucidated in the language of the 
many-body theory is the mechanism underlying a partial 
breaking of the quantized circulation in the presence of 
a rotating drive. 

We consider a system of N attractive bosons with mass 
771 on a ID ring with circumference 2'itR. Length and 
energy arc measured in units of R and h? / (2mR^) , re- 



spectively. The Hamiltonian for our system is given by 
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where ip^O) is the field operator, which annihilates an 
atom at position 9, and g (> 0) denotes the strength of 
attractive interaction. According to the GP mean-field 
approximation for the Hamiltonian (^, the ground state 
is either a uniform condensate or a broken-symmetry 
bright soliton, depending on whether the parameter gN 
is below or above the critical value, gN = 1. In contrast, 
all eigenstatcs of the original Hamiltonian are transla- 
tion invariant, and many-body theory predicts that the 
ground state is either a single {gN ^ 1) or fragmented 
{gN > 1) condensate |3|- 

Figure H (a) shows the low-lying spectrum obtained 
by exact diagonalization of the Hamiltonian ^. The 
dramatic change in the landscape of the energy spec- 
trum around gN ~ 1 is a consequence of the quantum 
phase transition between a single condensate and a frag- 
mented one. Figure n(b) presents the Bogoliubov spec- 
trum obtained from the Bogoliubov-de Gennes equations. 
By comparing Figs. 1(a) and (b), we find that the Bo- 
goliubov spectrum has a one-to-one correspondence with 
the many-body spectrum for gN < 1. For gN > 1, 
however, the many-body spectrum becomes much more 
intricate than the Bogoliubov one. In the Bogoliubov 
spectrum for gN > 1, there appears a Goldstone mode 
A' (the translation mode of the soliton) associated with 
the symmetry breaking of the ground state, the breath- 
ing mode B' , and the second harmonic of the breathing 
mode C". In Fig. Q] (a) for gN > 1, in contrast, a num- 
ber of quasi-degenerate levels appear with the density of 
states peaking around the Bogoliubov levels; we denote 
the corresponding groups as A, B, and C, respectively. 
The basis states for the diagonalization are restricted to 
the angular-momentum states I = 0,±1 (Zc = 1) un- 
less otherwise stated, and the field operator is given by 
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FIG. 1: (a) Excitation spectrum obtained by exact diago- 
nalization of Hamiltonian (1) for N = 200. The inset shows 
the corresponding result obtained with truncation Ic ~ 2 near 
the critical point, (b) Bogoliubov spectrum corresponding to 
(a), where branch A' represents the Goldstone mode, B' the 
breathing mode of a bright soliton, and C' the second har- 
monic of B' . (c) Energy gap AE between the ground and 
the first excited states in the many-body spectrum versus the 
total number of atoms A*' with gN — 1.4 held fixed. Triangles 
and circles denote results obtained with Ic = 1 and Ic ~ 2, 
respectively. 



i/'(6') = (co + e*^ci + e^'^c„i)/-/27r, where q is the an- 
nihilation operator of a boson with angular momentum 
hi. The validity of this cutoff has been confirmed by the 
inclusion of higher angular-momentum states as shown 
in the inset of Fig. ^ (a) with Zc = 2 |9(, where the energy 
landscape and the degree of degeneracy are unchanged 
from those oi l^ = 1- 

We denote the eigenstates of the Hamiltonian as \C)cr, 
where C is the total angular momentum. The index 
a = A,B,- ■ ■ labels the eigenstate in the ascending or- 
der of energy {E^c)^i < E^c)^ < ' ' ') for each C The 
states \C)a and \ — C)cr are degenerate in the absence 
of a rotating drive. At gN = 0, the eigenstates are 
the Fock states \C)a = \ni,n^i), where n±i denote the 
numbers of atoms with angular momenta ^ = ±1 and 
C = ni — n_i. The energies of the states |ni,n_i) 
are given by ni + n_i = j, and these states are thus 
(j -|- l)-fold degenerate. The index a = A,B,--- corre- 
sponds to the number of ^ = ±1 pairs being given by 
(j - |£|)/2 = 0, 1,---. For < gTV < 1, the energy 
branches are characterized by j, and the {j + l)-fold de- 



generacy is almost maintained. As gN approaches 1, 
each branch begins to ramify, and the energy landscape 
for gN > 1 is characterized by the index a and C as 
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There is no Goldstone 



mode because the ground state possesses the translation 
symmetry, and the lowest excited states | ± 1)^ have a 
finite energy gap AE = S|±i) ^ — ^|o>a' since the system 
is finite. However, the density of states above the ground 
state becomes higher for larger N, and the gap AE col- 
lapses as 1/N [Fig. n](c)]. The ground state is therefore 
unstable against excitations of the quasi-degenerate low- 
lying states. 

We construct the many-body counterparts of the 
bright soliton \^g) and the Goldstone mode \^e), such 
that (^el^e) = by superpositions of the ground and 
quasi-degenerate states: 



|*e) = e 
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(3o\0)A + J2l^ci\C)A + \-C)A] 
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<i>,)^-|vI/,H-^e-^'' ^ /:/3ai>C)A - I - /:)a) ,(3) 
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where L = J ddt/j^ {6){— {80)^(6) is the angular momen- 
tum operator, and Pcs satisfy J^c \Pc\'^ = 1- "^^"^ '^^" 
ergy cost associated with the superposition ^, -E'|*e> ~ 
^|ou = E£ l/3£l'(£^|£).4 - ^|o)J, isontheorderof 1/iV. 
This indicates that the symmetry breaking from the ex- 
act ground state \0)a to the bright soliton |^e) costs 
little energy, and the superposition thus occurs by an 
"infinitesimal" perturbation on the order of 1/A^. 

The emergence of quasi-degenerate levels also plays a 
crucial role in breaking the quantized circulation. In the 
presence of a rotating drive with angular frequency 2fl, 
the many-body ground state is either a single or frag- 
mented condensate, depending on whether /(g, A^, il) = 
{l-gN)/2-2{n-[n + l/2])'^ is positive or negative [lo| . 
where [D, -\- 1/2] denotes the maximum integer that does 
not exceed fl + 1/2. Figure |21 shows the total angular 
momentum of the ground state £g and low-lying eigen- 
values of the Hamiltonian H — 2nL + fl^ in the rotating 
frame. There clearly appear two distinct regimes where 
the density of states of excitations is low (/ > 0) and 
high (/ < 0) in the spectrum, and the circulation hCg/m 
is quantized only when / > 0. When the density of 
states above the ground state is sparse, the ground state 
cannot make a transition to states with higher angular 
momenta even if Q, is increased, since the term —2Q.L is 
not large enough to make up for the excitation energy 
due to a large energy gap. The total angular momentum 
£g therefore does not increase with Jl, and is quantized 
at integral multiples of N for / > 0. However, a number 
of branches with higher angular momenta decrease in en- 
ergy as / becomes negative; then, as fi is increased, one 
branch after another takes the place of the ground level 
upon intersection [Fig. [3 (b)]. The angular momentum 
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FIG. 2: Low-lying spectrum (thin lines, left scale) and ex- 
pectation value of the total angular momentum (thick line, 
right scale) of the ground state versus the angular frequency 



of the rotating drive for g ■ 
enlargement of (a). 



: 2 X 10"^ and N = 100. (b) is an 



of the ground state £g thus increases stepwise each time 
the excitation energy collapses. The interval between col- 
lapses of the excitation energy, i.e., the width of the steps 
of >Cg, becomes narrower as N becomes larger, eventually 
resulting in the breaking of the quantized circulation, as 
indicated by the slopes in the thick line in Fig.|3(a). The 
emergence of quasi-degenerate levels also induces symme- 
try breaking in a manner similar to Eq. |(2Jl. In fact, the 
ground state in the GP theory is a localized soliton for 

/<oIIi. 

We next investigate the effect of a symmetry-breaking 
perturbation by diagonalizing the Hamiltonian H + eV 
where H is given by Eq. iQJ and V = J dO'ip'^ (9) cos 6tl>{9) . 
As the perturbation is switched on, the ground state 
begins to localize by superposing the quasi-degenerate 
states. The largest eigenvalue X^ of the single-particle 
density matrix is plotted as a function of eN'^ in Fig. O 
The fragmented condensate approaches the single con- 
densate as eN^ increases. We note that AJ^ depends 
only on eiV^, which indicates that a single condensate 
is realized for e > N~'^. The fragmented condensate 
with iV 3> 1 is hence very fragile against the symmetry- 
breaking perturbation. 

The condensate fraction A}^ and the distribution of /?£ 
in Eq. can be derived analytically with the Bogoliubov 
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FIG. 3: The largest eigenvalue of the reduced single-particle 
density matrix obtained by the exact diagonalization (Aj^ ) 
and Bogoliubov theory (Ag ) versus e'N'^ for gN — 1.4. Inset: 
Distribution of \Pcf = U(/:!*(e)>P with e = 1 x lO"-*. The 
solid curves depict Eq. with Eq. (|^. 



approximation 1 1] . Since the frequency of the transla- 
tion mode of the broken-symmetry soliton is proportional 
to e^' ^ and small, the number of excited atoms in this 
mode is much larger than in the other modes. We thus 
consider only the translation mode, and obtain the deple- 



tion as (27V)-ie-i/2[(2i^) 



-1/2 _pl/2 



0{e)\ = 1-A 



(e) 

B ' 



where F = 7{3gN + A)^^^/[{5gN + 2){2gN - 2)^/^] is a 
decreasing function for gN > 1 . The behavior of Ag is in 
excellent agreement with the numerical result as shown 
in Fig. 13 Assuming that Pc\^)a in Eq- 12)1 corresponds 
to the projection of the Bogoliubov ground state on the 
angular momentum C, i.e., I3c\^)a ^ / |f e'(^"-^)*|*B), 
we obtain the distribution of Pc as 



\Pc 



1 



' V2^ 
,, 2(giV^l) 
7gN 



exp 



(2F) 



2d2 

1/2 , 



Fe 



1/2 



0(e) 



(4) 

(5) 



which is in excellent agreement with the numerical results 
as shown in the inset of Fig. |2| The behavior of the 
width d ex N^/^e^^^ indicates that the ccnter-of-mass 
fluctuation is proportional to N~^^'^e~'^^^. 

Finally, we investigate what happens to the exact 
ground state under the action of a quantum measure- 
ment jl2| . The ground state \0)a is prepared as an initial 
state with the number of atoms A^init so that gA^init > 1- 
Suppose that one atom is detected at position 6j in the 
j-th measurement. The postmeasurement state I^P^^^) is 
related to the premeasurement one |vl/(j-i)) by j^^-'^) = 

1|^{ej)\^^^-^^)/^{^U-^)\lP^0j)i;i0J)\^U'^)). The nor- 
malized single-particle density is given by n^^' {9) = 
(^(j)|^t(6»)^(6))|*(j))/7Vj, and the position of the (j + 1)- 
th measurement Oj^i is probabilistically determined after 
the j'-th measurement according to the probability distri- 
bution n'^^\0). Since each run is a stochastic process, we 
numerically perform sequential runs of independent sim- 
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FIG. 4: Ensemble-averaged largest eigenvalue Xj^ of the 
reduced single-particle density matrix versus the number of 
quantum measurements j for gNinit = 1.6. The solid curve 
denotes the analytical result Q. Inset: Ensemble-averaged 
value of \l3cf = U{/:|*^-''')P after j = 5,20 and 50-time 
measurements. The solid curves depict Eq. ||1J with (f = 
4i(<7iVi„it - l)/(7.giVi„it). 



ulations and take the ensemble average for each j. We 
find that the ensemble-averaged value of the condensate 
fraction AJ^ is independent of TVinit for a fixed .giVinit, 
and that it monotonically increases as shown in Fig. 0] 
Therefore, when iVinit ^ 1, AjJ reaches the order of one 
for j/iVinit <C 1, indicating that the fragmented conden- 
sate rapidly becomes a single condensate by the action 
of quantum measurements. The single-particle density 
rS^'>{6) localizes and the state |^g ) reaches the form of 
Eq. ©. 

Suppose that the initial state is a uniform superposi- 
tion of the sohton state / d9oA{9Q)(j)'^^[9 — Oq), where 6*0 
is the center of mass and A{9o) is a constant for j ~ 0. We 
can show [ll| that the probability distribution is changed 
by the j'-time measurements as A^^Oq) oc exp(— 4cj/0q) 
with c = 2(giVinit — l)/(7giVinit), and hence the distri- 
bution of |/3£p reduces to the Gaussian form |0J with 
(P = 2cj, i.e., d ex -y/j for j < VMnit- The centcr-of- 
mass fiuctuation thus reduces to A^o oc j"^^^ by the 
j'-time measurements. The condensate fraction can also 
be derived analytically [ll| as 
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1 



[1 



c + ce 



-2a 



+V(1 - c + ce-2")2 - 4c(l - 2c)(l - e-")2 ,(6) [g] 



where a = l/(8cj). These results are in excellent agree- 
ment with the numerical ones as shown in Fig. 01 

In conclusion, we employed the exact diagonalization 
method to investigate the simultaneous emergence of 
a bright soliton and the Goldstone mode in a ID at- 
tractive BEG. We found that the existence of a num- 
ber of quasi-degenerate states above a critical strength 
of attractive interaction makes the ground state fragile 



against symmetry-breaking perturbations. By introduc- 
ing a symmetry-breaking perturbation or quantum mea- 
surements, we showed that a localized state with an en- 
hanced condensate fraction is generated by a superposi- 
tion of quasi-degenerate many-body states, and that the 
resulting state closely resembles the Bogoliubov ground 
state. The mechanism underlying a partial breaking of 
the quantized circulation is also elucidated in terms of 
the many-body spectrum. 
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